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We explore very fine scales of scalar dissipation in turbulent mixing, below Kolmogorov 
and around Batchelor scales, by performing direct numerical simulations at much finer 
grid resolution than is usually adopted in the past. We consider the resolution in terms 
of a local, fluctuating Batchelor scale and study the effects on the tails of the probability 
density function and multifractal properties of the scalar dissipation. The origin and 
importance of these very fine-scale fluctuations are discussed. One conclusion is that 
they are unlikely to be related to the most intense dissipation events. 



1. Introduction 

Direct numerical simulations (DNS) of the equations governing turbulent phenomena 
have now become an important tool in both research and applications. Faithful results 
require the equations to be solved on a grid of adequate resolution. Most solutions of 
the past adopt a standard resolution based on an average measure of the smallest scale 
estimated from dimensional considerations. Our purpose here is to show, for the case of 
turbulent mixing, that such standard measures do not resolve a range of very fine scales 
that are present in reality. We demonstrate the essentials of this feature by performing 
several numerical simulations whose resolution is much finer than that adopted routinely, 
and comment on the practical relevance of the scales missed in past calculations. 

We focus on the practically important case of Sc = vjn > 1, where Sc is the Schmidt 
number, v the kinematic viscosity of the fluid and k the scalar diffusivity. The smallest 
scale of mixing is generated by the balance of turbulent stretching and diffusion (Batch- 
elor 1959), and is related to the Kolmogorov scale rj = (v z /e) 1 ^ 4 through the identity 
rj B = rj/\/~Sc; here, e is the average value of the instantaneous (and local) energy dissi- 
pation rate e of the turbulent kinetic energy. The expression for e reads as 



where Uj(x, t) is the velocity fluctuation in the direction i. In simulations of a fixed size, 
in which the advection diffusion equation and the Navier-Stokes equations are simulta- 
neously solved, the largest scale has traditionally been maximized by resolving no more 
than rj B . Resolving this scale is also the goal of most experimental efforts, though they 
often fall short for large Peclet numbers Pe = u'L/k, where v! and L, respectively, are 
the characteristic large-scale velocity and length scale of the flow. 

However, since the time it was realized that the local energy dissipation rate e has a 
multifractal character (e.g., Sreenivasan & Meneveau 1988), and that its average value 
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Figure 1. (a) Energy spectrum in the far-dissipation range for runs with N = 128 and N = 512, 
as indicated in the legend. In (b) we have fitted the data by E(k) ~ k a cxp(— ck) where k = krj, 
with a = 3.3 and c = 7.1 as in Chen et al. (1993). 

e is therefore not adequate for quantifying small-scale characteristics (Grant, Stewart & 
Moilliet 1962), it has been clear that, locally, length scales smaller than rj appear in a flow. 
This is obvious from the definition of the local Kolmogorov scale r/(x, t) — (i^ 3 /e(x, i)) 1 / 4 , 
in which e(x, t) is a highly fluctuating variable: the larger the value of e locally, the 
smaller is the local Kolmogorov scale 77, and vice versa. In view of the relation between 
the local values of the Batchelor and Kolmogorov scales, namely 77s (x, rj) = 77 (x, t) / \/~Sc, 
it becomes similarly apparent that the existence of very small values of n may lead to t\b 
values that are much smaller than the average value rj B . This intuitive (and incomplete) 
reasoning will be supplemented by further explanation momentarily, but it is enough to 
note here that such fine scales have not been explored before. 

Section 2 presents a brief description of the simulations and resolution effects, and the 
main results are stated in section 3. This is followed by a discussion in section 4 of how 
the finest dissipation filaments are generated. Conclusions are stated in section 5. 

2. Resolution criteria and the far-dissipation range 

The advection-diffusion equation for the passive scalar is solved together with the 
Navier-Stokes equations for the homogeneous and isotropic turbulent velocity field. The 
simulations are carried out within periodic boxes of fixed size L 3 = (27r) 3 . The flow 
is maintained stationary by a random volume forcing. The passive scalar is rendered 
stationary through a mean scalar gradient in the y-direction. Except for enhanced reso- 
lution, about which some comments follow, the pseudo-spectral methods used here are 
quite standard (Schumacher & Sreenivasan 2003; Vedula, Yeung & Fox 2001). 

For DNS of box-type turbulence on an iV 3 grid based on pseudo-spectral methods 
control of alias errors requires that modes with wavenumber higher than k m = V2N/3 
be truncated (Patterson & Orszag 1971). Specifically, this removes double and triple 
aliases in three dimensions. The usual resolution criterion is expressed as a number for 
k m rj which, by extension to cases with Sc > 1, is 



A proper resolution of the smallest scale requires that the global minimum of tjb be 
resolved by the grid which can be written as 



kmVB — 1-5 • 



(2.1) 



min X;t [77B(x,/j)] 
A 



> 1 . 



(2.2) 
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Table 1. Parameters of the numerical simulations. A is the (equidistant) grid spacing. The 
integral length scale for the velocity field is L u = 7r/(2u§) J °° dkE(k)/k and that for the scalar 
field is Lg = n/(28 2 ) J °° AkEe(k)/k; here E(k) and Eg(k) are the spectral densities of the 
velocity and scalar fields, respectively, u% is the mean square of the velocity in the x-direction, 
and 6' 2 that of the passive scalar. T av /TE is the averaging time in units of the large scale eddy 
turnover time Te = 3u x 2 /2e. 



By comparing (2.2) with the usual resolution criterion l|2.1|) . we get 

n r>— ~0.5, (2.3) 

where we used A = L/N = 2tt/N and k m = \/2N/3. It is apparent that the criterion 
H2.2(l is stronger. Our simulation will satisfy (|2.2|l and its results will be compared with 
data from those with nominal grid resolution. 

Unfortunately, the demands of this fine resolution restrict us to rather low values of 
the Reynolds number (R\ < 24, see table 1 for details). This necessitates that most 
of our flow scales are in the viscous range of turbulence. Figure ^a) shows the energy 
spectrum for R\ = 10, evaluated with two resolutions, N — 128 and 512. It is seen 
that the noise floor is reached some 30 orders of magnitude below the peak signal. In 
figure ^b), we verify that the energy spectrum in the far-dissipation range falls off as 
E(k) ~ k a exp(— cfc) with k = krj. For the lower Reynolds number case, one gets a = 3.3 
and c = 7.1 as did Chen et al. (1993) for the same range of dissipative wavenumbers; this 
range is indicated by the two vertical dotted lines in the figure. 

3. Resolution effects on the scalar dissipation field 

Figure 2 shows that the small scales of the scalar dissipation field eg , defined by 

appear as filamented structures in a planar cut. This is especially highlighted in the 
three-dimensional rendering of figure [3] which shows that very intense parts of scalar 
dissipation rate appear as sheets. Box counting of such events confirms a dimension close 
to 2. While no immediately discernible differences are apparent between the two scalar 
fields obtained with conventional and the present ultra-fine resolution (see the two upper 
panels of figure 2), the differences in eg can be detected more easily at several positions. 
This is seen quantitatively in figure 01 which shows that the resolution matters far away 
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Figure 2. A slice in the (x-z) plane of the scalar field (upper row) and the corresponding scalar 
dissipation field (lower row in logarithmic units, black for maximum and white for minimum) 
at Sc — 32 and Rx = 10. Left column: low resolution case with N = 128. Right column: high 
resolution case with N = 512. A quarter of the plane is shown with an area of 13Vq B x 13l77 s . 



from the mean, or for the tails of the probability density function (PDF) of eg. For 
eg <C e$, the PDF scales with an exponent of 1/2 and the curves for the two resolutions 
do not coincide. Significant differences in the large-amplitude tails (eg /eg > 20) are also 
apparent. An analytical result for the tails of p(eg) exists for Pe — > oo in a smooth white- 
in-time flow. Using the Lagrangian approach, Chertkov, Falkovich & Kolokolov (1998) 

1 /3 

and Gamba & Kolokolov (1999) deduced the behaviour to be p(eg) ~ exp(— e g ). Our 
finding here seems to be consistent with this result though our data are in the Eulerian 
frame. 

Another result that emphasizes the presence of very fine scales is shown in the left 
panel of figure which plots, from three well-resolved simulations, the PDF of the local 
Batchelor scale ??b(x) — rj(x)/^fSc generated by mixing. The data for all three values 
of Sc are advected in exactly the same flow. rj B is indicated for each Sc by a vertical 
dashed line. It is clear that scales substantially smaller than rj B do exist, and that they are 
inaccessible to the standard resolution l|2.1|l . The right panel of figure |3] illustrates that 
an increase of the Reynolds number will broaden the range of »7s(x) and thus magnifying 
the effect. 

To highlight the sensitivity to fine-scale resolution, we have calculated the generalized 
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Figure 3. Isosurfaces of the scalar dissipation field for R\ — 24 and Sc = 32. The level is 
z = eg /eg = 11. Isosurfaces form flat and curved sheets which was also confirmed by a box 
counting analysis of isolevel sets. 
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Figure 4. PDF p(z) of the scalar dissipation field z = eg /eg for A*' = 128 and N = 512 at 
Sc — 32. The outer panel shows the PDF in log-log scales while the inset shows a log-linear 
plot. The predicted tail behavior for a smooth white-in-time flow in the limit of Pe — > do 
(Chertkov al. 1998) is indicated in the inset. For the latter fit we included only the data with 
z = eg /eg > 20, as indicated by the dotted vertical line in the inset. 



dimensions 

Dq{q) = lim J_ kgX>?(r) ( j 

where r is the diameter of the subvolumes Bi(r) of the successive coarse graining and 
the measure Hi(r) = e^B~(^j/ee (Hentschel & Procaccia 1983). In the mid panel of fig- 
urc|21we compare the generalized dimensions for two runs, taking the same scaling range. 
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Figure 5. Left: Probability density function (PDF) of local fluctuations of the Batchelor scale 
in the mixing problem. The Taylor microscale Reynolds number is 10 and the Schmidt numbers 
are 2, 8, and 32. The computational domain is 512A on the side, where A is the grid spacing. 
The vertical dashed lines close to the maximum of each PDF indicates the average Batchelor 
scale rj B . The dotted vertical line corresponds to one A. Right: Reynolds number dependence 
of the fluctuations of the local dissipation scale, ??(x). The Batchelor scale for Sc = 32 and the 
grid spacing A are indicated by dotted lines. Since all lengths are rescaled by rj, both these 
parameters collapse for the two Reynolds numbers. 



Differences for q < are readily apparent. This part of the D q curve is dominated by 
low magnitudes of scalar dissipation. Clear differences are visible in the spatial distri- 
bution of regions of scalar dissipation below a certain small threshold value for the two 
runs (compare left and right panels). D q (q) for q > also shows differences, with the 
inadequately resolved data slightly underestimating the peak dissipation regions. Larger 
effects of inadequate resolution correspond to low amplitude regions rather than to high 
amplitude regions. They can be identified with wavcnumbers k > k* where k* is the 
wavenumber at which the scalar dissipation spectrum eg(k) — 2nk 2 Eg{k) peaks. 

It is not surprising that poor resolution — which, in some sense, translates to increased 
noise — has a stronger effect on regions of low eg than those of high eg where the signal- 
to-noise ratio is effectively high. 



4. Scales of very high and very low scalar dissipation 

We note that two different situations of large Pe are of interest: (a) Sc = 0(1) at high 
Reynolds number and (b) Sc^> 1 but at lower Reynolds numbers. In laboratory exper- 
iments dealing with liquid flows, the Reynolds numbers are moderate and the Schmidt 
numbers quite high (Sreenivasan 1991; Buch & Dahm 1996; Villermaux & Innocenti 1999; 
Catrakis et al. 2002). Most air experiments (Sreenivasan 1991a; Mydlarski & Warhaft 
1998) are at higher Reynolds number and Sc = 0(1). Numerical simulations have consid- 
ered either moderately high Re and modest Sc (Vedula, Yeung & Fox 2001; Watanabe & 
Gotoh 2004) or low Re and large Sc (Yeung, Xu & Sreenivasan 2002; Brcthouwer, Hunt 
& Nieuwstadt 2003; Schumacher & Sreenivasan 2003; Yeung et al. 2004). Our situation is 
neither (a) nor (b) exactly, and we may expect small-scale scalar fluctuations influenced 
by the forcing which arises from velocity fluctuations. 

The fluctuations of the velocity field around the Kolmogorov scale — in the intermediate 
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Figure 6. Generalized dimensions D q of the scalar dissipation field for nominal and ultra- fine 
resolutions. The values of D q , shown on the left panel, were obtained by fits in the range 
r £ [L/2 6 , L/2 4 ] with L = 2tt. The right panels show (without any additional numerical smooth- 
ing) slices through a snapshot for the data set eg/~eg < 0.025 illustrating the differences in the 
calm regions that contribute to D q for q < 0. 



dissipation range — have been discussed in Frisch & Vergassola (1991). A set of velocity 
increments, having a certain Holder exponent h in the inertial range — i.e. 5^v = |u(x + 
£) — u(x)| ~ £ h with 8tv being the velocity increment over a generic scale £ in the inertial 
range — and occupying a spatial support with a (fractal) dimension D(h), is associated 
with a dissipation length n(h). The "roughest" increments scale with h m in & n d cause 
the largest spikes of energy dissipation persisting down to the smallest dissipative scale, 
V(hmin)- The extent of the intermediate dissipation range is basically the width of the 
probability density functions as plotted in figure^] (right). A global minimum of the local 
Batchelor scale would be given by the condition 

minMx,i)] = %^, (4.1) 
x .* VSc 



which means that strongest scalar dissipation is controlled by T](h m i n ) = f 3 / 4 / max(e) . 
To see if such a scale is related to the most intense scalar dissipation events, we estimate 
ee ~ k(9 2 ) / £ 2 by setting £ = n(h m i n ) / 'y/Sc which is smaller than rj B and obtain 



max(e e ) ~ 6 2 \l — ■^- ± ■ (4.2) 

This connects directly the regions of high scalar dissipation to locations of large energy 
dissipation rate, as anticipated in section 1. Limitations of this expectation can be seen 
in figure |7| which plots the joint PDF, p(e, eg). Large values of e are not necessarily 
connected to large values of eg. The crossing point of the dotted lines which is marked 
as a black square indicates that an extreme event as given by l|4.2l) is not present. 

A second possibility is that the most intense scalar dissipation events are at scales in the 
viscous-convective range. Taking scalar increments SgQ over this range, i.e. rj B < £ < rj, 
one gets due to Batchelor (1959) the expression 

v 



S e 6 2 = (0(x + £)- 0(x)) 2 ~ 1~g\l — log ^— ) . (4.3) 
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equal increments of the logarithm to base 10, decreasing from 0.5 in steps of —0.5. The crossing 
point of the dotted lines is the maximum of eg following 14.21 . The dashed line indicates the 
maximum if that maximum were to occur in the viscous-convective range. 



Rx = 10 R x = 24 

A fc = a k ± ib k z < 10~ 3 z > 42 z < 10~ 3 z > 69 z > 42 z > 11 

(01,02,03) 23% 87% 28% 90% 82% 58% 

(o±i6,-2o) 77% 13% 72% 10% 18% 42% 

Table 2. Eigenvalue analysis of the velocity gradient tensor diUj at sites where z = eg /eg 
exceeds/falls below a threshold. The mean fraction of each of the possible eigenvalue solutions 
is given for two Reynolds numbers and a Schmidt number of 32. Ai + A2 + A3 = 0. 



A scale-dependent maximum of 

ee(l) pb k— ~nee\j- ~ 2 (4.4) 

can be calculated via dee(£)/d£ = 0. This occurs at I* = rj B y/e, which is larger than 
rj B . The corresponding value eg(l*), shown by the dashed line, remains below the dotted 
horizontal line given by Ij4.2|l . The related scale is larger than the Batchelor scale which 
eases somewhat the strong resolution requirement 1)2.2(1 . 

Of particular interest is the nature of the velocity field near the most intense or the 
least intense parts of the scalar dissipation. This was examined through the eigenvalue 
analysis of the velocity gradient at sites where z = eg /eg exceeded or fell below a chosen 
threshold. Trends with Sc were discussed by Schumacher & Sreenivasan (2003), but here 
we focus only on the highest Schmidt number. The results of table 2 show that the 
pure straining motion (three real eigenvalues) and rotation (complex conjugated pair of 
eigenvalues) contribute to low amplitude scalar dissipation. It may be thought that the 
scale I of maximum dissipation is related closely to the most compressive strain events, i.e. 
I = y/ k/ max(|7|). The present data for all levels of z, some of which are reported in table 
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2, indicate that this is not so. Indeed, a broad range of compressive rates — not merely 
the maximal magnitudes — are associated with intense scalar dissipation. This is already 
apparent from the analysis of Ashurst et al. (1987) (see their figure 9b which reports 
an isotropic flow for Sc=0.5). While the most compressive eigenvector was preferentially 
aligned with the direction scalar gradient, corresponding events of maximum eg were 
preferentially aligned at an angle of about 20° . We explicitly note that this is the Eulerian 
point of view, and that a Lagrangian analysis following incipient fronts to the stage of 
their maturity may yield a different result. 



5. Discussion 

The possibility that the resolution requirement could be more stringent than is conven- 
tionally believed has been discussed to some detail in Sreenivasan (2004), but the details 
outlined in the present paper have not been explored before. When there is a significant 
overlap of the intermediate dissipation with the viscous-convective range, extreme values 
of the scalar dissipation are determined by the roughest velocity increments of the iner- 
tial range of turbulence. An important example in which these fine scales would make 
a difference is non-premixed turbulent combustion (Bilger 2004). There, the scalar dis- 
sipation rate of the mixture fraction enters as a basic quantity, e.g. for the modelling of 
jet-diffusion flames. Chemical reactions take place at the stochiometric mixture fraction 
in sheets of sub-Kolmogorov thickness. One expects steep gradients across such layers 
and strongly varying scalar dissipation. These variations are not captured in the flamelet 
equations where only the statistical mean enters the expansion parameter (Peters 2000). 

A broad example where resolution effects can be important is the multifractal scaling 
of e and eg. The spiky structures in space and time, which are most prominent in the 
dissipative scales, are thought to affect appropriate turbulent field even for scales larger 
than rj or vb, as appropriate (Frisch 1995). 

To summarize, the very fine scalar filaments that were resolved here do not seem to 
be associated with the most intense scalar dissipation. The small-scale stirring in the 
flow seems to interrupt a further steepening process which is known as the formation 
of mature fronts. Similar findings were made for two-dimensional turblence at Sc = 1 
(Celani et al. 2001). Clearly, the Reynolds number of the advecting flow has an impact 
on this issue simply because the intermediate dissipation range might "overshadow" the 
viscous-convective range completely for sufficiently large R\. A more conclusive study of 
this issue will require higher Reynolds numbers and will be part of future work. 

Computations were carried out using the NPACI resources provided by the San Diego 
Supercomputer Center and on the Jiilich Multiprocessor (JUMP) IBM cluster at the 
John von Neumann- Institute for Computing, Jiilich. Special thanks go to Herwig Zilken 
(Jiilich) for his help with figure|21 PKY and KRS acknowledge support by the US National 
Science Foundation and JS from the Deutsche Forschungsgemeinschaft. 
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